Tests for non-randomness in quantum jumps 
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In a fundamental test of quantum mechanics, we have observed 228 000 quantum jumps of a 
single trapped and laser cooled 88 Sr + ion. This represents a statistical increase of two orders of 
magnitude over previous similar analyses of quantum jumps. Compared to other searches for non- 
randomness in quantum mechanical processes, using quantum jumps simplifies the interpretation of 
data by eliminated multi-particle effects and providing near-unit detection efficiency of transitions. 
We measure the fractional reduction in the entropy of information to be < 6.5 10 -4 when the 
value of any interval between quantum jumps is known. We also find that the number of runs of 
successively increasing or decreasing interval times agrees with the theoretically expected values. 
Furthermore, we analyze 238 000 quantum jumps from two simultaneously confined ions and find 
that the number of apparently coincidental transitions is as expected. Finally, we observe 8400 
spontaneous decays of two simultaneously trapped ions and find that the number of apparently 
coincidental decays from the metastable state agrees with the expected value. We find no evidence 
for short- or long-term correlations in the intervals of the quantum jumps or in the decay of the 
quantum states, in agreement with quantum theory. 

PACS numbers: 03.65.Ta 42.50.Lc 



An axiom of quantum mechanics is that we cannot 
predict the result of any single measurement of an ob- 
servable of a quantum mechanical system in a superpo- 
sition of eigenstates. Testing this principle is important, 
not only for basic science, but also for applications such 
as quantum random number generators (QRNG's) and 
potential quantum computers. It is therefore surpris- 
ing that in spite of the many experiments sensitive to 
quantum mechanical effects, only a few experiments have 
explicitly searched for non-random behavior in long se- 
quences of repeated quantum measurements. In 0,0, 
the randomness of the path of a single photon after a 
beamsplitter was used to build QRNG's. In Q, the ar- 
rival times of decay products of unstable nuclei were used 
to test the statistics of quantum decay. Although both 
these systems rapidly give excellent statistics, detector 
inefficiencies limit the conclusions that can be drawn re- 
garding the unpredictability of quantum mechanical mea- 
surements. Furthermore, both systems are insensitive to 
certain types of non-random behavior: averaging over 
many particles in a collection of nuclei could obscure non- 
random behavior of single systems, patterns in emission 
times or photon arrival times could be overlooked be- 
cause of inefficient detectors, and because beamsplitters 
are always somewhat biased, QRNG's based on them are 
designed to be insensitive to consecutive runs of transmis- 
sions or reflections. All these problems can be avoided by 
observing the times of quantum jumps in a single atom [4j 
because transitions between atomic levels can be detected 
with near-unit efficiency with no multi-particle effects [j| . 
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FIG. 1: Partial energy level diagram, transitions and lifetimes 
in 88 Sr + . We drive the 422-nm transition to Doppler cool and 
detect the ions, and the 1092-nm transition to prevent optical 
pumping into the D 3 / 2 state. A laser of bandwidth < 10 kHz 
drives the 674-nm transition to produce quantum jumps to 
and from D 5/ / 2 states. 



Here, we analyze the intervals between quantum jumps 
of a trapped ion and answer the question "after we have 
measured the length of one interval, how much have we 
decreased the uncertainty of the value of any interval we 
subsequently measure?" This question is directly relevant 
to quantum logic operations on any quantum mechanical 
system |9| ; if the result of one logic operation affects that 
of subsequent (and supposedly independent) operations 
then the operation of the gate is compromised. It is also 
critical for demonstrating the suitability of a quantum 
mechanical system as a random number generator. 

Earlier work |6j examined 640 quantum jumps in a 
single 198 Hg + ion confined in a Paul trap |7|. Ref. [8j 
later reported a different but limited analysis of 10 000 
quantum jumps. In a substantial increase in statistics 
over this previous work, we analyze 228 000 quantum 
jumps from a single ion comprising continuous data sets 
of ~10 000 events each. We take advantage of our greater 
statistics to also test for unexpected correlations be- 
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FIG. 2: Number of 422-nm photons counted by the detector 
over time, when all of the transitions shown in Fig. are 
driven simultaneously. 

tween transition times of multiple ions, analyzing 238 000 
quantum jumps and 8400 spontaneous decays of pairs of 
trapped ions. 

We confine 88 Sr + ions in a linear Paul trap an d si- 
multaneously drive the transitions shown in Fig. Fig- 
ure |2 shows a sample of the measured 422-nm scattering 
rate as a function of time. It displays the well-known 
characteristics of quantum jumps: the ion rapidly scat- 
ters resonant 422-nm light when it is in the S1/2 <-> 
P1/2 D3/2 manifold, but not at all when it is in one 
of the metastable D 5 / 2 magnetic sub-levels. The scatter- 
ing rate changes abruptly whenever a 674-nm photon is 
absorbed or emitted. According to quantum theory, the 
exact times of these changes, and the intervals between 
them, should be unpredictable. 

To test this assumption, we analyze the number of 
422-nm photons that scatter from the ion and reach a 
photomultiplier tube in a measurement time t meas (typi- 
cally a few ms, followed by a 200 /is dead time) when all 
three transition in Fig.^are simultaneously driven. Each 
data set of approximately 10 000 measurements is taken 
over approximately 30 minutes. From these data we ob- 
tain four series U = {u\, U2, 113, . . . , itzv} of sequential 
time intervals, defined as follows. We label as "bright" 
the series of intervals during which the ion continuously 
scatters 422-nm light, and "dark" the series of intervals 
during which the ion fluorescence is continuously absent. 
In addition we analyze the series of intervals between 
successive emissions ("emission") and absorptions ("ab- 
sorption") of a 674-nm photon. 

Before each data set, we adjust the intensities of the 
laser beams to control the rate of quantum jumps. The 
bright intervals are distributed exponentially, with aver- 
age values ranging from 70 to 200 ms, while the aver- 
age values for dark intervals range from 30 to 150 msec. 
When integrated over the long times of a data set, the 
distributions of dark intervals deviate slightly from a 
purely exponential form for interval durations greater 
than 100 msec. This is due to drifts in the frequen- 
cies and intensities of laser light at the site of the ion, 
which cause slight variations in the average bright and 
dark interval durations. In addition, in the presence of 
a degeneracy-breaking magnetic held B (we vary B from 



10~ 6 T to 2 10~ 4 T), transitions from each D 5 / 2 Zeeman 
sub-level contribute distinct exponential distribution to 
the distribution of dark intervals. In spite of their dif- 
ferent responses to these two effects, the dark and bright 
intervals in our experiment give statistically identical re- 
sults. We conclude from this that these effects neither 
introduce nor mask correlations between intervals on the 
time scales of our statistical tests. 

Previous studies of randomness in quantum mechanics 
have used several algorithms for answering the question 
regarding predictability of intervals posed at this begin- 
ning of this paper: for example, analyzing the distribu- 
tion of sums and differences of intervals 6] , verifying that 
the probability of particle decay is constant over short 
and long time-scales |3J , and analyzing lengths of contin- 
uously increasing and decreasing runs of intervals 0, |(| . 
We calculate the difference between the entropy and con- 
ditional entropy of sets of intervals, because this directly 
answers our question, and because it results in an up- 
per limit on the association of adjacent and non-adjacent 
pairs of intervals. 

To do this, we first divide our data sets into continuous 
strings of 1000 events and normalize each interval time by 
the average interval time of its string. This nearly elim- 
inates effects due to slow drifts of the rate of quantum 
jumps. Using each interval Ui only once, for k > 1 we tab- 
ulate the number of occurrences N m ^ n of pairs {itj, Ui+fe} 
for which m € [t m ,t m+ i) and u i+k £ [t n ,t n+ i). We 
choose t m f n \ at intervals that result in roughly uniform 
values of N mn ^> 1 for all {m, n}, although we find 
that our results do not depend on the degree of unifor- 
mity. We set the number of bins in each row and column 

ip^max — Tlmctx) to be 13 . 

Next, following Press ^2] we ca U Pm- the probability 
that the first interval in a pair (m) lies within the range 
[t m ,t m +i); an d p n the probability that the second in- 
terval (ttj+fc) lies in the range [t„,t n +i)- We also call 
Pm.n = N m , n /N (where N is the total number of pairs) 
the probability that the pair falls in bin {m, n}. The 
entropies are then defined as 

H(x) = -^p.„ In p. n (1) 

n 

m 

and the conditional entropies are then defined as 

H(x\y) = -y]p m ,n ^ (Pm,n/Pm-) (3) 

H(y\x) = - ^Pm.n In {Pm,n/P-n) ■ (4) 
m,n 

These define two joint measures of the association be- 
tween the first and second intervals of any pair: 

H(x)-H{x\y) ... 
U W y) = HM (5) 
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U(y\x) = 



H(y) - H(y\x) 

H(y) 



(6) 



U(x\y) represents the amount of information that the 
distribution of first intervals gives about the distribution 
of second intervals an d vice versa for U(y\x). For 
an infinitely large data set, if the intervals in a pair are 
statistically independent of each other then U(x\y) = 
U(y\x) = 0, and if they are completely associated with 
each other, U(x\y) = U(y\x) = 1. However, statistical 
fluctuations in the N^s limit how closely U(x\y) and 
U(y\x) approach 0. 

Figure |2| shows U(x\y) and U(y\x) as functions of the 
number of intervals k that separate each pair {ui,Ui + k} 
for the four different interval types. The values do not 
depend on k. Figure0]shows the average values of U(x\y) 
and U(y\x). They are consistent with each other and 
with the value 0.00062 that is calculated if all of the N m ^ n 
were randomly distributed about the average value q = 
N rn^ ax with a distribution width yfq. These facts show 
that intervals are uncorrelated when they are separated 
by 1 to 20 intervals, which corresponds to approximately 
2 ms to several seconds. 

As a complimentary test for short- and long-term or- 
der, we search for deviations from the expected distri- 
bution of consecutively increasing or decreasing runs of 
elements of length m. As an example, the subsequence 
{. . . , 3, 15, 11, 5, 7, 2 . . .} contains a run down of three el- 
ements (15, 11 and 5) and a run up of two elements (5 
and 7). If the data set U has a total of N elements, the 
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FIG. 3: The joint functions U(x\y) (solid lines) and U(y\x) 
(dashed lines) from Eqs. E|and|S|as a function of the number 
of intervals between pairs for bright and dark intervals and 
intervals between emissions and absorptions of 674-nm pho- 
tons. Heavy dashed lines indicate the average values for each 
interval type. 
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FIG. 4: The average values of U (x\y) (solid boxes) and U (y\x) 
(hollow boxes) from Fig. [3] 



expected number of runs up or down of length £ is 



n (e 2 + 1 - i) + £ (4 - e) + 1 

(J+2)! 



(7) 



This equation assumes that no element in U is repeated. 
However, we measure the intervals Uj in integer multiples 
of t meas . Often a run is terminated by a repeated inter- 
val value, although it is equally likely that the run should 
continue. We account for such cases by tallying two run 
lengths, each with half the weight of a single run. We 
make similar adjustments when an interval occurs three 
times in a row. We combine all of the data to obtain 
N runs (£) and their deviation from the expected values, 
shown in Fig. EI Because the N runs (£) are not indepen- 
dent of each other, a \ 2 test does not apply. However, 
each point in the figure is within reasonable agreement 
with the expected values. 

A critical assumption in quantum algorithms involving 
multiple quantum systems is that in the absence of inten- 
tionally applied coupling, the individual systems do not 
communicate with each other. We do not expect to ob- 
serve cooperative effects in the transitions times of pairs 
of atoms. However, other workers have seen many more 
apparently simultaneous quantum jumps 13] and spon- 
taneous decays [lj] than expected. We take advantage 
of our statistics and search for signs that transitions in 
ions stored simultaneously in the trap are correlated. 

We analyze 248 000 quantum jumps in two ions sep- 
arated by > 20 /xm, and count the number of times in 
which both ions appear to change state during the same 
measurement time t meas . In such events, the number of 
detected 422-nm photons 71422 in one measurement time 
is less than a threshold consistent with no ions fluoresc- 
ing, and is immediately followed or preceded by a value 
of 71422 that is greater than a threshold consistent with 
two ions fluorescing. The probability of such events is de- 
termined by the probability per unit time for a single ion 
to change states, the average 422-nm photon scattering 
rate per ion, the two threshold values and the measure- 
ment time t meas . We also account for the possibility of 
misinterpreting the scattering rate from a single ion for 
that of two ions due to insufficient resolution between the 
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FIG. 5: Differences between the expected and measured val- 
ues of the number of runs up (squares) and down (triangles) 
for the four interval types. The values of the number of runs 
of length t > 6 have been magnified, and the runs up and 
down have been slightly offset along the x-axis for clarity. In 
total, each data set has ~ 57 000 runs of various lengths both 
up and down. 

count rate distributions of one and two ions. The total 
number of coincidental transitions into the D 5 / 2 states is 
expected to be 308, and we measure 320. Also, a total 
of 316 coincidental transitions out of the D 5 / 2 states is 
expected from our data, and we measure 313. In addi- 
tion, we find that the observed numbers agree with those 
produced in Monte Carlo simulations of the data. 

We also analyze the spontaneous decay of simultane- 
ously trapped ions, which may be more sensitive to ef- 
fects such as interactions with background gas. A brief 



(< 0.2 s) saturating pulse of 674-nm light excites the 
atoms to the D 5 / 2 states while the 422-nm light is absent. 
After the 674-nm light pulse, the 422-nm light is returned 
and we monitor the 422-nm photon scattering rate every 
tmeas = 5 msec. We observe 8400 decay processes that 
start with two ions in the D 5 / 2 state and finish with no 
ions in the D 5 / 2 state. In total, 26 of these transitions ap- 
pear to occur during the same measurement time t meas . 
From a measured decay rate of 410 msec in our system, 
we expect to see 24 (4) of these processes. This, too, is 
consistent with the behavior of the ions being random, 
and agrees with the results of [Till 

In conclusion, while it is impossible to prove random- 
ness, we have seen no signs of non-random behavior over 
short and long time scales after analyzing 228 000 quan- 
tum jumps in single ions, 238 000 quantum jumps in two 
simultaneously trapped ions, and 8400 spontaneous de- 
cays of two ions. Practical QRNG's and quantum com- 
puters would use fewer quantum interactions than those 
analyzed here. The present sensitivity is sufficient to 
show that these applications are not affected by corre- 
lations due to non-randomness of quantum mechanics. 
However, at times short compared to the transit time of 
light across the ion 0] , and times long compared to the 
decay rate of the atomic state 01 ; t ne Weisskopf-Wigner 
approximation breaks down so that the distribution of 
decay times is no longer exponential. In this case, we 
would expect to see non-random behavior in the interval 
times of quantum jumps. Observations in these regimes 
are out of the reach of modern experiments, but are an 
intriguing possibility for future work. 

This work was funded by DOE through the LDRD 
program. We would like to thank Malcolm Boshier for 
carefully reading this manuscript and for valuable dis- 
cussions, and Richard Hughes for initially bringing this 
topic to our attention. 
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